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Abstract 

We revisit the evaluation of one-loop modular integrals in string theory, employing new 
methods that, unlike the traditional 'orbit method', keep T-duality manifest throughout. In 
particular, we apply the Rankin-Selberg-Zagier approach to cases where the integrand func- 
tion grows at most polynomially in the IR. Furthermore, we introduce new techniques in 
the case where 'unphysical tachyons' contribute to the one-loop couplings. These methods 
can be viewed as a modular invariant version of dimensional regularisation. As an example, 
we treat one-loop BPS-saturated couplings involving the d-dimensional Narain lattice and 
the invariant Klein /-function, and relate them to (shifted) constrained Epstein Zeta series of 
0(d, d; Z). In particular, we recover the well-known results for d = 2 in a few easy steps. 
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1 Introduction 

In essence, closed string theory is a quantum field theory of infinitely many fields, obtained 
by tensoring two infinite towers of left-moving and right-moving excitations subject to a 
level-matching constraint |TJ. As a result, any closed scattering amplitude at one-loop is an 
integral over two parameters, the Lagrange multiplier T\ G [— I, h] for the level-matching 
constraint and the Schwinger time T2 > parameterising the loop. Due to diffeomorphism 
invariance on the string world-sheet, the identification of the proper time is not unique and 
the integrand F(t) is a modular function of the complex parameter t = T\ + 1T2 (in general 
not holomorphic). To avoid an infinite over-counting, the domain of integration is restricted 
to a fundamental domain T for the modular group T = SL(2,Z), thereby removing the 
ultraviolet divergences from the region T2 — > which usually arise in quantum field theory. 
The usual field theoretical infrared divergences from the region T2 — > 00 are in general still 
present, due the existence of massless particles in the spectrum]*] 

In general, computing such 'one-loop modular integrals' is a daunting task, (in part) due 
to the unwieldy shape of the domain T . For specific amplitudes describing BPS-saturated 
couplings in the low-energy effective action however, the integrand simplifies and the mod- 
ular integral can be computed explicitly. One of the simplest instances occurs for certain 

1 In the following, we restrict our attention to closed string theories without physical tachyons, but allow for 
'unphysical tachyons', i.e. relevant operators which do not satisfy the level-matching condition, such as those 
present in heterotic models. 
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anomaly-related couplings in the ten-dimensional heterotic string theory \2\. In this case the 
integrand is the elliptic genus F = <P(t), a weak holomorphic modular form with a singu- 
larity at the boundary of T , and the modular integral can be evaluated by applying Stokes' 
theorem j2}[3). 

In lower dimensions, however, the low-energy couplings depend non-trivially on the 
geometric moduli of an internal d-dimensional torus T d . Indeed the integrand function is 
typically of the form ^(^JXj+fc^/ where 3>(t) is again a weak holomorphic modular form 
and r^ d+k/d ^ is the partition function of the Narain lattice associated to the torus compacti- 
fication. The integrand is not holomorphic, so the integral cannot be reduced to a line in- 
tegral on the boundary of T via Stokes' theorem. Rather, the main technique for dealing 
with such modular integrals in the physics literature has been the 'unfolding trick' or 'orbit 
method', pioneered in (4}|5| and generalised in many subsequent works p|[6}|T4). In a nut- 
shell, this method consists in extending the integration domain J 7 to a simpler region (the 
strip n G [-1/2, 1/2], t 2 > 0, or the full upper half plane X\ e R,t 2 > 0) at the cost of 
restricting the sum over momenta and windings to suitable orbits of lattice vectors. While 
leading to a very useful expansion at large volume, this method has the drawback of ob- 
scuring the invariance of the resulting low-energy coupling under the automorphism group 
0(d + k, d;Z) of the Narain lattice. Although in some simple cases it is still possible to 
rewrite the result in terms of known automorphic forms of 0(d + k,d;Z), this is in general 
not easy to achieve. For k = 0, = 1, it was conjectured in p5| that the result of the one- 
loop modular integral could in fact be expressed (in several different ways) as a constrained 
Epstein zeta series, manifestly invariant under 0(d, d; Z). 

The purpose of this note is to introduce a procedure for evaluating modular integrals 
which keeps manifest any (additional) symmetry of the integrand function. To this end we 
apply the Rankin-Selberg-Zagier (RSZ) method, a close relative of the 'unfolding method' 
which has been a standard tool in the mathematics literature (see e.g. |16] for a survey). In 
short, the main idea is to insert a non-holomorphic Eisenstein series in the modular integral 
of interest,^] 




gcd(m,n)=l 



compute the latter 'deformed' integral by applying the 'unfolding trick' to the sum over m, n 
for large 3?(s), and finally obtain the desired integral by analytically continuing the result 
to s = (where the Eisenstein series is known to reduce to a constant). Due to the growth 
of F(t) as T — > ioo in cases of physical interest, the integral in is infrared divergent, 
while the original Rankin-Selberg method was restricted to functions F(r) of rapid decay at 
the cusp. Fortunately, the Rankin-Selberg method was extended by Zagier in |17) to allow 
for functions of moderate growth at the cusp, by introducing a hard infrared cut-off T 2 < T 
on the Schwinger time, unfolding the sum over m, n at finite T, and giving a prescription 

throughout the paper, d}i = r^ 2 drj dT2 denotes the SL(2;Z) invariant measure on the hyperbolic plane, 
normalised so that fjrdfi = n/3. 
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for renormalising the integral as the infrared cut-off is removed. With this renormalisation 
prescription, one can view the replacement as a stringy analogue of dimensional reg- 
ularisation, which preserves modular invariance^] Moreover, the integral for s ^ 
literally arises for certain BPS couplings, e.g. the D 4 R 4 couplings studied in p0| . Using the 
RSZ method, we shall evaluate the renormalised modular integral fj- djiF^^ exactly for 
any value of d, in terms of a constrained Epstein zeta series of 0(d,d,Z), thereby proving 
the conjecture in [15|. We shall also recover the celebrated result of [6| for d = 2 in just in a 
few easy steps, illustrating the power of this approach. 

While the Rankin-Selberg-Zagier method outlined above is very efficient for functions 
F(t) of polynomial growth, which is typically the case for BPS couplings in type II string the- 
ory, it is unfortunately inadequate for functions F(t) of exponential growth, which typically 
arise in heterotic amplitudes, due to the ubiquitous unphysical tachyon. More specifically, 
we are interested in modular integrals of the form 



dv&(T)r {d+k/d)/ (i.2) 



where ^(t) is a weak holomorphic modular form of weight w = —k/2 with an essential 
singularity at the cusp, <£(t) ~ e -2niKT _|_ £)(i) with k > (for heterotic strings, k = 1 but 
our method works equally well for any k > 0). Since the Ti -average of the polar part of 
F(t) = ^(t) r^+jy) vanishes, the divergence of the integral ( |1.2| is not worse for k > than 
for k = 0, but the RSZ method nevertheless fails and one must resort to different techniques. 

In the second part of this work we shall develop a new procedure for dealing with the 
above class of modular integrals, which relies on representing the holomorphic part ^(t) of 
the integrand in terms of a Poincare serie^j 

*(t)= E <K7-t), (1.3) 



for a suitable function i/?(t), and then applying the unfolding trick to the sum in ( |1.3| . Our 
procedure is close in spirit to the Rankin-Selberg method, in particular it keeps manifest the 
0(d + k, d;Z) symmetry of the Narain lattice, however note that it does not involve any 



auxiliary modular function for unfolding the fundamental domain as in (1.1 1, rather it uses 
the function F itself or part of it. 

A naive implementation of this idea, however, is hampered by the fact that the Poincare 
series of a modular form of non-positive weight is not absolutely convergent and thus its un- 
folding is not justified. We circumvent this problem by deforming ^(t) to a non-holomorphic 



3 Other modular invariant infrared regulators have been proposed, e.g. Jl8) . It is also possible to regulate the 
integral by subtracting the non-decaying part of F(t), as in (61 . As we shall see, it is straightforward to relate 
these different regularisation schemes. For an early use of hard infrared cut-off regularisation methods in string 
theory, see e.g. (T9[ . 



Here, Too C T is the stabiliser of the cusp ioo, generated by the triangular matrices ? = (g \) ,n e 
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Poincare series E(t\, t 2 ; s) such that 



«J>(T)=limE(Ti,T 2 ;s), (1.4) 

s— s-0 

applying the unfolding trick for large 3?(s) and recovering the desired integral in the limit 
s — > 0. To illustrate our method, we shall compute the modular integral J djiFu^ j(r) and 
represent it in terms of a 'shifted constrained Epstein zeta series', which (unlike the treatment 
in |3j[8j) makes its invariance property under 0(d, d,Z) manifest. 

Before closing this introduction, we note that the RSZ method has already been useful in 
string theory for studying the distribution of the graded degrees of freedom in tachyon-free 
oriented closed string vacua, and their connection to the one-loop free energy [21 1. Using 
this method it was shown that non-tachyonic string configurations are characterised by a 
spectrum of physical excitations that not only must enjoy asymptotic supersymmetry but 
actually, at very large mass, bosonic and fermionic states are bound to follow a universal 
oscillating pattern, whose frequencies are related to the non-trivial zeroes of the Riemann £- 
function. Similar studies have then been generalised to higher genus 1 22 - 24j where similar 



constraints on the interactions of physical states are expected to emerge. 

Furthermore, we anticipate that further development of the techniques presented in this 
work will allow to compute a variety of other modular integrals of interest in string the- 
ory. In particular, in vacua with broken supersymmetry, it would interesting to study the 
behaviour of the one-loop effective potential near points of symmetry enhancement, where 
extra massless states appear. These typically lead to singularities in the vacuum energy that 
are difficult to analyse using the standard 'unfolding method' 1 25 [26 [, whereas they should 



be fully captured by our new approach. This would allow one to probe the stringy behaviour 
around the points of enhanced gauge symmetry. It would also be interesting to generalise 
our methods to higher genus. 

The outline of the paper is as follows. In Section |2j we recall the original Rankin-Selberg 
method for functions of rapid decay, and its generalisation to functions of moderate growth 
by Zagier. Our paraphrasing of |l7| is mainly due to the 'basic identity' which appears in 
the body of the proof of [17] which we later require. In Section|3j we apply the RSZ method 
to one-loop modular integrals of symmetric lattice partition functions and recover the cele- 
brated result of [6j in very few steps. We further give a proof of a conjecture in 1 15 [, clarifying 



the relation between the constrained Epstein zeta series of |T5| and the Langlands-Eisenstein 
series studied in |20[ . In Section |4j we develop a variation on the method of Rankin-Selberg 
and Zagier for functions of rapid growth at the cusp but with finite, or at most power-like 
divergent, modular integrals. We apply our procedure to the modular integral of symmetric 
lattice partition functions times the modular /-invariant (and its images under the action of 
the Hecke operators), and we express it in terms of a novel shifted, constrained Epstein Zeta 
series of 0(d, d; Z). The appendix collects some properties of Kloosterman sums used in the 
text. 
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2 A brief review of the Rankin-Selberg-Zagier method 
2.1 Rankin-Selberg method for functions of rapid decay 

Assume that F(r) is an automorphic function of the complex variable r = X\ + iT2 of rapid 
decay at the cusp, i.e. such that F(t) vanishes faster that any power of t 2 at t 2 -¥ 00. The 
Rankin-Selberg transform TZ*(F,s) of F is defined as the Petersson product 

7T(F,s) = \?{2s) E / dp j FfaTjO, (2.1) 
M)=i 

between F and the non-holomorphic Eisenstein series 

E*(T; S )^r(2s)£(r; S ) / E(t;s) = \ £ g = £ [3 ( 7 • t)] s . (2.2) 

M)=i 

Here £*(s) = 7r _s,/2 F(s/2) £(s) is the completed zeta function with simple poles at s = 1,0 
and zeroes in the strip < 3?(s) < 1. The integral over the fundamental domain T can be 
unfolded on the strip S = {t 2 > 0, — \ <X\ < \ }, so that 



7r(F;s)=r(2s) / **g*F(r)=?(2s) f° dr 2 t 2 s - 2 F (r 2 ) , (2.3) 

XT s JO 



F (t 2 ) = / 1/2 dnF(T). (2.4) 

J-1/2 



dri dr 2 

3 

where F (t 2 ) is the constant term of F, 

-1/2 
1/2 



Thus, 7Z*(F, s) is proportional to the Mellin transform of Fq. Now, £*(t; s) is well-known to 
be a meromorphic function in s, with simple poles at s = 0, 1, satisfying the first Kronecker 
limit formula 

£ * (T; S) = 2(^1) + K 7 " log(47r T2 l?/(T) |4) ) + ° {S ~ 1} ' (2 - 5) 

where 7 is the Euler-Mascheroni constant and t](x) = c^'^YXn^iX ~ <?") is the Dedekind 
function (as usual, q = e 2mT ). These statements follow from the Chowla-Selberg formula 

£*(t;s) = C(2s)t{ + T(2s- 1)^+2 £ |NrW 2s(N ) r 2 1/2 K s 1 (2tt|N|t 2 ) e 2 ^ , 

N^O 2 

(2.6) 

where (N) is the divisor function 

a t (N)= £ d f , (2.7) 
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and Kt (z) is the modified Bessel function of the second kind. The properties 



r(s) = r(l-s), K t (x) =K- t (x), <T t (n) = n f <r_ t (n) , 
ensure that £*(t;s) satisfies the functional equation 
E*(t;s) = E*(t;1-s). 



(2.8) 



(2.9) 



It then follows that lZ*(F;s) inherits the same analytic and functional properties of E*(t;s), 
z'.e. it is a meromorphic function of s with simple poles at s = 0, 1 and symmetric with respect 
to the critical axis K(s) = \, 



7T(F;s) = n*(F;l-s). 



(2.10) 



This method was used in the mathematics literature 1 27 281 to establish analytic properties 



of certain Dirichlet L-series (see e.g. [ 16] for a survey). More importantly for our purposes, 
the residue at s = 1 of the Rankin-Selberg transform is equal to (half) the average value of F 
on the fundamental domain T , 



Res K*(F;s)\ a=1 = \J^d}iF= -Res ft*(F;s)| g=0 . (2.11) 
This in principle provides a way to evaluate the integral of the automorphic function F on 



the fundamental domain, from the Mellin transform of the constant term Fq [21-231. In 



particular, the residue at s = 1 depends only on the behaviour of Fo( T 2) near T2 = . 

2.2 Rankin-Selberg method for functions of moderate growth 

In physics applications, F is rarely of rapid decay. Fortunately, the Rankin-Selberg method 
has been adapted to the case of automorphic functions of moderate growth by Zagier in fll7| , 
which we paraphrase below. 

Let F(t) be an automorphic function whose behaviour at the cusp x = ioo is of the form 



F(t)~<p(t 2 )+0(t 2 - n ) (VN>0), 



where 



(2.12) 



(2.13) 



for suitable c, G C, a, G C, n,- G N. For this class of functions, following Zagier 1 17 [, we define 
the Rankin-Selberg transform as 

7T (F;s) = ?(2s) / dr 2 z^ 2 (F - <p) , (2.14) 

where Fq(t2) is the Ti-constant term ( 2.4 1 of F. The integral ( 2.14| | converges absolutely when 
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5R(s) is large enough (namely, 3R(s) > for all z). As we shall see^J lZ*(F;s) can be 

meromorphically continued to all s, with possible poles at s = 0, 1, ocj and 1 — a,, and is 
invariant under s i— > 1 — s. Moreover, (half) the residue of 1l*(F;s) at s = 1 gives a prescrip- 
tion of the (otherwise divergent) renormalised integral of F on the fundamental domain. 
To establish this, we shall use a combination of 'hard infrared cut-off and 'zeta function 
regularisation', i.e. consider the (manifestly finite) integral 

K T (F;s) = f d}iF{r)E k (r;s) (2.15) 
J Ft 

on the "cut-off fundamental domain" Tj = T H {t 2 < T}. It is a fundamental domain for 
the "cut-off Poincare upper half plane" 

H T = nn{T 2 <T} - |J s a/c , (2.16) 

(fl / c)ez 2 , 

c>l,(fl,c)=l 

where S a / C is the disk of radius 1/ (2c 2 T)) tangent to the real axis at a/c. Defining xt to be 
the characteristic function of Hj and performing the same unfolding trick as in ( |2.3| with 
F ■ xt i n place of F, we obtain 

7^(f;s) = T(2s) / drxd^F^)^- 2 , (2.17) 



Using ( |2.16[ ), this may be rewritten as 
( T 

f dr 2 rdnFtr)^- 2 - £ / dT!dT 2 F(T)T 2 s - 2 . (2.18) 

J0 J ~2 C>l,{a,c) = l, JS «l<: J 



K T (F;s) = C(2s) 



V 



a mod c 



Now, the disc S a / C is mapped to H n {t 2 > T} by any element 7 = | ] E T. For fixed 



a/c, all such elements are related by right multiplication by Foo. Thus, the last term in the 
bracket in ( 2.18| can be rewritten as 



/ d^F(T) £ [S( 7 -T)] s d^. (2.19) 

The sum over 7 reproduces the Eisenstein series E(t;s), modulo the term t| due to the 
restriction c > 1. Putting everything together, we find 

K T (F;s) = r(2s) / dT 2 F (T 2 ) t^ 2 - / d^F (E*(t;s) - £*(2s) t|) . (2.20) 



5 This statement can be seen right away by noticing that <p(T2) is annihilated by the differential operator 
□ = n'=i[A — 0ii{di — 1)]"' +1 , where A is the Laplacian on H, and applying the standard Rankin-Selberg method 
to the rapidly decaying function OF. We are grateful to D. Zagier for pointing this out. 
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The symmetry s i— > 1 — s may be restored by further subtracting the second constant term 
from E*(t;s), so that, after multiplying by 7r~ s r(s) and rearranging terms, 



/ dpF(Ti,T 2 )E*(T;s)dp+ / d^F(T 1 ,T 2 ) (e*(t;s)-E*(°)(t; ; 
^ T " r (2.21) 

/*/ /"CO 7 

= T(2s) / dT 2 F (T 2 )T 2 s - 2 -r(2s-l) / dT.FolT,)^ 1 - 5 , 
JO " JT 



where E*(°) (t; s) = £*(2s) t| + £*(2s — 1) T 2 s is the constant term in the Fourier expansion 
of E*(t;s). Since this is X\ -independent, the product F£*'°'(t;s) appearing in the first line 
can be replaced by FoE*(°)(t;s) without changing the result of the integral over the strip 
J 7 — Tj. Now, the terms in the second line evaluate to 



(■ I poo 

J dT 2 Fo(T 2 )T 2 s - 2 = 7r(F;s)/n2s)- J dr 2 (F - q>) r^ 2 + h T (s) , 

/CO f-CO 
dr 2 F (t 2 ) t^ 1 - 5 = J dr 2 (F - <p) t^ 1 - 5 + fc^(s) , 



(2.22) 



/T JT 
where /zj- and are incomplete Mellin transforms of q>, 

rT /*oo 

h T (s)= dT 2 <p(T 2 )r 2 s - 2 , ti T (s)= dr 2 <p(r 2 ) t 2 - x - s . (2.23) 
■JO JT 

A key fact about the class of functions q> in ( 2.13| is that their (complete) Mellin transform 
vanishes, therefore h'j{s) = —hj{l — s). Moreover, integrating q> t 2 ~ 2 once, 

l r . m (i\ni-m n-s+Ki- 1 i m -y- 



Using this and rearranging terms, we arrive at Zagier's basic identity, Eq. (27) in [17|, 

n*(F;s)=[ duFE*(r;s) + [ du ( F E* (r; s) - <p E*<® (r; s)) 

h T h-T T V > (2.25) 

-C(2s)Ms)-r(2s-l)Ml-s)- 

Evidently, the r.h.s. of ( 2.25| l is independent of T, meromorphic in s, invariant under 
s i — y 1 — s, and analytic away from s = 0, 1 (where E*(r;s) has a simple pole) and from 
s = a,, 1 — «; (where 7z7-(l — s), respectively hj{s), has a pole of degree n, + 1)0 Thus, when 
no a, coincides with 0, 1, 

6 The apparent pole at s = 1/2 cancels between the last two terms in l |2.25) . 
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where 3>(s) is an entire function of s. Moreover, the residue at s = 1 is 



Res ft*(F;s)| s=1 =-Res [£*(2s) h r (s)] s=1 - Res [£*(2s - 1) fe r (l " *)Ui 
+ | / d^F(ri / T 2 )+/ dp (F(t 1/ t 2 ) - <p(r 2 )) 

Letting <p(i" 2 ) be an anti-derivative of cpixj), 

a, i=-\ 



(2.27) 



(2.28) 



a,=l 



and defining the renormalised integral as 



R.N. / d^P(Ti / T 2 ) = / dpF(Ti,T 2 ) + / dp (F(T a ,T 2 ) - <p(T 2 )) - 0(7") 



lim 

T -S-oo 



/ d^F(T!,T 2 ) - f(T) 



(2.29) 



which is by construction T-independent, eq. ( 2.27} may be rewritten as 



R.N. J dp F(r lr T 2 ) = 2 Res [7T (F; s) + £*(2s) ft r (s) + C(2s - 1) fc r (l - s)] 8=1 - <p(T) . 

(2.30) 



In fact, the r.h.s. of (2.251 is itself the renormalised integral 



ft*(F;s) =R.N. ^d^F(T 1/ T 2 )£*(r;s) / 



(2.31) 



with F(t 1/ t 2 )E*(s / t) / <p(r 2 ) E*(°)(t;s) and £*(2s) ft r (s) + £*(2s - 1) M 1 ~ s ) playing the 
role of F{x\, t 2 ), <p(r 2 ) and <p(T) in ( |2.29[ ) r respectively. 

For functions of rapid decay the renormalised integral reduces to the usual integral and 
h-r,f{T) vanish, hence ( |2.30[ ) reduces to ( 2.11| . More generally if 3?(a,) < 1 for all i, the 
integral J^duF still converges, and one can take the limit T — > oo in ( |2.30} and recover ( 2.11} . 
If however one of the 9fJ(a,-) > 1, the integral Jj- Fdp diverges like <p(T) as a function of 
the infrared cut-off, and ( |2.30| l provides a renormalisation prescription which depends only 
on the divergent terms with 9fJ(«;) > 1 in ( |2.13| |. 

Of course, the renormalisation prescription ( 2.29} is by no means the only possible one. 
For example, one may decide to subtract the non-decaying part from F and integrate the 
remainder on the fundamental domain, defining 



R.N/ J dfiFfaTz) = ^dp (F(n,T 2 ) - <p([r 2 })) . 



(2.32) 



Here we denoted by [t 2 ] the imaginary part of 7 • t, where 7 is an element of T which maps T 
into the standard fundamental domain T (so [r 2 ] = t 2 if r £ J 7 ). The renormalised integrals 
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( 2.32 | and ( 2.29| l differ by a finite quantity 



4 = R.N. / d^F(Ti / T 2 )-R.N. / / dpiF{r lr T 2 ) = lim / dji (p([r 2 )) - f{T) 

which can be computed explicitly. For example, if all the n/s are zero, 

x r l\ !-«< . 3. 1\ 



, (2.33) 



^ = E< 



i=l 



CC; — 1 



(2.34) 



with d> ai b being the Kronecker symbol, and 2Fi(a,b;c;z) the standard hypergeometric func- 
tion. Notice that the case a, = 1 is well encoded in ( 2.34[ > since, in the limit a, — > 1, 



in 



1 1 - a, 3 1 
2' 2 ' 2' 4 



1-log 



3V3 X 



( a! -l) + 0(( a! -l) 2 ), 



(2.35) 



and thus A = 1 - log 



3 Lattice modular integrals and constrained Epstein zeta series 



In this section, we apply the Rankin-Selberg method to the evaluation of the integral of the 
lattice partition function^] 



A 11 E <' 



— 7TT2 .M g27riTj m; n' 



(3.1) 



on the fundamental domain. The left-handed and right-handed momenta 

Pl(r) i = ( mj ± (gij T Bi/) ) , 



(3.2) 



depend on the geometric data of the compactification torus, i.e. the metric gu of the T d and 
the NS-NS two-form Bu, that together parameterise the symmetric space 0(d,d)/0(d) x 
O(d), also known as the Narain moduli space. In the last equality of ( |3.1| , 



M 2 = { mi + B^g^nij + B fl n l ) + n^n' = vl + v\ 



(3.3) 



is the mass-squared of a string ground state with Kaluza-Klein momentum m, and winding 
number n l , with being the inverse metric. The lattice partition function is manifestly in- 
variant under 0(d, d; Z), but also under SL(2; Z) T . This invariance is exposed after Poisson 
resummation with respect to m,-, 



T (d,dM' B ) = V det g E ex P 

(m' ,n')i 



(m 1 + n i T)g ij {mi + n)T) t # 

-7r \-2mBijtn w 

T 2 7 



• (3-4) 



7 Here and in the following we set a! = 1, and we suppress the explicit dependence of the lattice on X\ and Tj. 
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The partition function for the Narain lattice clearly belongs to the class of functions consid- 



ered in Section 2.2 since, in the limit T 2 — > 00, 



~ Tr 



d/2 



(3.5) 



that matches eq. ( 2.13| with with I = 1, a, = d/2 and ft,- = 0. Using Zagier's extension of 
the Rankin-Selberg method, as summarised in the previous section, we can then compute 
the renormalised integral 



K*(r {d4 y,s) =R.N. J^d F r {d4) (g / B)E*(s / T) / 



(3.6) 



and especially its residue at s = 0, which is proportional to the IR finite one-loop integral of 
the lattice partition function, 



I d (g,B) =R.N. J^r m {g,B) 



(3.7) 



Before proceeding with the explicit calculation of the integral ( |3.6| , we notice that the 
Rankin-Selberg transform TZ*(r/ did ^,s) is an eigenfunction of the Laplace operator acting on 
the moduli space 0(d, d) /0(d) x O(d) 



A SO (W) K*{r m ,s) = l(2s-d)(2s + d-2)n*(r m ,s), 



with 



I a ! 7 



d gi, 



(3.8) 



(3.9) 



This follows straightforwardly from the following differential equations satisfied by the par- 
tition function of the Narain lattice |l5| and by the Eisenstein series 







A som - 2 A SL(2) + \ d(d - 2)1 r m (g t B) , 



(3.10a) 



A SL (2)-2s(s-l) E*(t;s), 



(3.10b) 



where 



A 



1 2 f d 2 a 2 

SL(2)-2T 2 ^ + 3^ 



(3.11) 



is the Laplace operator on the hyperbolic plane. 
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3.1 Constrained Epstein zeta series in dimension d > 2 

For a generic d-dimensional lattice ( p.l} , 



M<0 = ^ r , ^(T 2 )= T2 } (3-12) 

s + 4 d - 1 log t 2 if a = 2 , 



2 

and as expected, the integral 



l^dpr m (g r B)E*(s r r) (3.13) 

is power-divergent, logarithmically divergent and absolutely convergent for s + \d — 1 > 0, 
s + j d — 1 = and s + j d — 1 < 0, respectively. In all cases, however, the renormalised 
integral \3.6\ is finite and is given by the Mellin transform 



Tl*(r m ;s)=C(2s) J dT2T s+d/2-2 e -nr 2 M e 2nir im -n 

" rrii,n' 

C*(2s)r(s + U-i) . , 



(3.14) 



Here, m • n = m, n ; , and we have denoted by £ v {g, B; s) the constrained Epstein zeta series 
in the vectorial representation of 0(d, d; Z), introduced in |l5| 



4(«,b; S )-E'%#, (3-15) 

m,n J 1 

which converges absolutely for s > d (as usual, a primed sum does not involve the contri- 
bution from nii = n' = 0). It is useful to define the completed constrained Epstein zeta series 



£ d v * (g, B; s) = n~ s r(s) r (2s - d + 2) (g, B;s) , 
so that 

7i*(r (w) ;s) = 4*(g / B;s + Jd-l). 
From eq. ( 2.10| | it follows that £ v *(g, B;s) satisfies the functional equation 

£i*(g,B;d-l-s), 



(3.16) 



(3.17) 



(3.18) 



in agreement with eq. ( |3.8[ ). 

These properties, together with the invariance of £ v *(g,B;s) under the ring of 0(d, d)- 
invariant differential operators [15|, implies that the constrained Epstein zeta series coin- 
cides with the degenerate Langlands-Eisenstein series for 0(d, d) based on the parabolic 
subgroup P with Levi subgroup R + x SO(d — l,d — 1) poj 29 30). Moreover, from the 
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general statement below ( 2.25 1, it follows that, for d > 2, £y(g, B;s) has simple poles at 
s = 0,d/2 — l,d/2 and 1, as indicated in [29|. 



Using the functional equation ( 3.18} , we see that ( 3.14} is equivalent to 
R.N. ^ dfi r m (g,B) E*(s, t) = 4* (g, B;\d-s) . 
For d > 2 one can easily extract the residue at s = 1 to get 



(3.19) 



(3.20) 



where in writing the last expression we have made use of the functional equation ( |3.18| . 

Notice that the first equality in ( 3.20| establishes Theorem 4 in |T5| rigorously. Moreover, 
comparing with Eq. (C.2) in |20] (and dropping the superfluous volume subtraction), we 
recognise that the constrained Epstein zeta series is in fact equal to 



jSO{d4) ■ 



(3.21) 

where E^^ M J is the Langlands-Eisenstein series, introduced in [20 . This relation can also 
be checked by comparing the large volume expansions given inEq. (C.7) in |20) and in 
Appendix C.l of p}. 



3.2 Low dimension 

The cases d < 2 are special, since the integrals are at most logarithmically divergent, and the 
delta-function constraints in the definitions of the constrained Epstein zeta functions can be 
explicitly solved. 

For d = 1, the lattice partition function reduces to 



r (1/1) (R) = ^ £ e -7rT 2 [( m /R)2+(, lR )2] e 2inx l mn i 



(3.22) 



with <p(T2) = y/T2. The modular integral l\ is finite and coincides with the renormalised one. 
The constrained Epstein zeta series ( |3.15| evaluates to 



4(£,B; S )=2C(2s) (R 2s + R- 2s ), 
and thus 

fc*(r (u) ; 8 ) = ty(g,B;s - |) = 2? (2s) T(2s - 1) (r 1 ^ + R 2 *" 1 



(3.23) 



(3.24) 



The r.h.s. has simple poles at s = 0, 1 and a double pole at s = 1 12, in agreement with the 
general statement below ( |2.25| |. The residue at s = 1 produces then the standard result for 
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the modular integral 
h 



n 



d F r (u) (R) = ^ (r + r- 1 



(3.25) 



that is invariant under R 1-4 1/R as required by the 0(1,1;Z) symmetry of the lattice. 
Clearly the same result arises by unfolding the sum over m, n in ( |3.22| >. 

For d = 2, the modular integral f-p dji £(2,2) i s logarithmically divergent as T — > 00. The 
Rankin-Selberg transform is however finite for large K(s) and still given by the constrained 
Epstein zeta series (3.151. Once more, the constraint can be explicitly solved to get the stan- 
dard expression for the integral. To this end, it is convenient to parameterise the two-torus 
in terms of the complex structure modulus U = U\ + HI2 and Kahler modulus T = T\ + 1T2, 
so that 



PL 



m 1 + Um 2 + T(n 2 -Un 1 ) 



PR 



nil + Uni2 + T(n 2 — Un 1 
V2T2IT2 



(3.26) 



Now, we notice that the most general solution of the constraint m\n x + w^n 2 = is given by 
the elements in the disjoint union 



(m 1 ,m 2 ,n 1 ,n 2 ) e <Si U5 2 , 
with Si, S2 being the sets: 
«5i = {(m 1 ,m 2 ,0,0) , (m 1 ,m 2 ) E Z 2 } , 

S2 = {(cthi,cm2, —drh2,drh\) , [ih\,m%) e Z 2 ,^pA.{ih\,m%) = 1 , (c,d) e Z ,d > 1} . 
The contribution of the first set to £y*(T, U; s) easily gives 

T 2 U 2 



(3.27) 



(3.28) 



C(2s) E 

(mi,m2)GZ 2 
(mi,m)/(0,0) 



I mi + Uni2\ 2 



21lC{2s)E*(U;s), 



(3.29) 



where £*(t; s) is the SL(2, Z) invariant Eisenstein series (2.2 1. For solutions belonging to the 
second set, the mass-squared factorises as 



\Pl\ 



\PR\ 2 = 



\rhi + Urh2\ 2 \c + Td\ 



U 2 



so that the contribution of £2 gives 

T 2 



r(2s) £ 



\c + Td\ 2 



E 



T 2 



U 2 



(3.30) 



\ih\ + Ufh2\ 



2E*(U;s) E 



T 2 



Tdl 2 



rf>i 



(3.31) 
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These two results combine into a simple expression for the constrained Epstein zeta series, 
n*{r {2i2 y,s) = £l*(T,U-,s)=2E k (T-,s)E k (U-,s). (3.32) 

This relation is actually a consequence of the group isomorphism 

0(2,2;Z) ~ SL(2;Z) T x SL(2;Z) U x Z 2 , (3.33) 

and of the decomposition of the vectorial representation of 0(2, 2) in terms of the bi-fundamental 
(2,2) of SL(2) x SL(2). 

It is clear from ( |3.32| l that £y*{T, U,s) has a double pole at s = and s = 1. Once more, 
this is in agreement with ( |2.25| l, since, upon using £* (s) = 1 / (s — 1) + | (7 — log 47r) + O (s — 
1) and h(T) = T s /s, the second line in this equation has a double pole at s = 0, 1, 

r(2s) h T ( S ) + r (2s - 1) mi - = - 2 { S -\y + 2(^1) [1 °s (47rT) - 7] + • • • • (3 - 34) 

To compute the renormalised one-loop integral 1%, one then needs to extract the residue at 
s = 1. Using the first Kronecker limit formula (|2.5|l, one finds 



ft*(%2);s) 



1 



+ 



2(s-l) 2 s-1 L 



7-ilogfl67t z T 2 Lr 2 |)7(T)i7(Lr)| 4 y + ... . (3.35) 



Combining ( 3.34| , ( |3.35[ ) with ( |2.30| l and using <p(7~) = log(T), one arrives at the follow- 
ing expression for the renormalised integral of the two-dimensional Narain lattice partition 
function 



l z = R.N. l T r m (T, U)dfi = - log (471 e-T T 2 U 2 |//(T) ^(U)| 4 ) . 



(3.36) 



We stress that Eq. ( |3.36 l is the result of the renormalisation prescription ( 2.29| l. It differs by 
a finite constant from the renormalised integral computed in |6f using the renormalisation 
prescription (232 1. Indeed, using ( |2.34| one arrives at 



\ T (r m (T, U) - t 2 ) d H = - log (^4^T 2 U 2 |i/(T) 7(11) 



3^ 



(3.37) 



in agreement with J6|. In our opinion, the derivation of 1 3.37 1 via the Rankin-Selberg-Zagier 
method is considerably simpler than the original derivation in (6J. 

Finally, we comment on the properties of the Rankin-Selberg transform ( 3.32| l and the 
renormalised integral ( |3.36[ ) of the lattice partition function 1% 2 ) under the action of T- 
duality invariant differential operators. For d = 2, Eq. ( 3.10a| l becomes 



-2A (T) +A (T) +A (U) 



r (2,2) = / 



(3.38) 



(T II) 

where A^'J are the analogue of the Laplacian ( |3.11[ ) acting on the hyperbolic T, U-plane. 
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However, one may check that -T(2,2) satisfies the stronger equations 



a£L r, 



A 



(T) 



A 



(U) 



SL (2) 1 (2,2) - ^ S L(2) 1 (2,2) - ^ S L(2) 1 (2,2) • 



(3.39) 



Combining these equations with ( |3.10b[ ) r we find that TV* (1(2,2)} s ) must be an eigenmode of 

a sl(2) anc ^ a sl(2) se P ara tely, with the same eigenvalue \s{s — 1). This is of course manifest 
from the explicit result ( |3.32[ ). Moreover, due to the subtraction of the second order pole 
in (335}, the renormalised integral is a quasi-harmonic function of the moduli, namely it 
satisfies 



AT) r _ AU) j 
A SL(2) *2 ~ A SL(2) J 2 



(3.40) 



This is again manifest from the explicit result ( |3.36| |. These considerations will become fruit- 



ful when computing integrals with unphysical tachyons in Section 4.4 



3.3 Decompactification 

It is useful to ask about the behaviour of the modular integral l\3.6) when the radius R of one 
circle in the cf-dimensional torus is sent to infinity. For simplicity, we restrict to the subspace 
of the Narain moduli space where the lattice partition function factorises into 



B ) = r (i,i)( R ) x r (d-l,d-i)(g/B)/ 



(3.41) 



where g, B are the metric and NS-NS two-form on the remaining (d — 1) -dimensional torus. 
To investigate the limit R — > oo, we consider the modular integral with a hard cut-off, and 
unfold the lattice sum r^iy Following the same reasoning as in eqs. ( 2.18| |-( Z20| | one finds, 
for large enough K(s), 



n* T (r m ;s)=R I dji J"((j-i /( i-i) {g, B; t) E*(t;s) 



J-l/2 t£ Vm/Q 
00 rV2 d Tl dT 2 



T J-l/2 x\ 



_ 7tR 2 |ib-ht| 2 

Re T 2 | r, 



r (d-i,d-i)(g' B;t)E*(t;s) 

\d-i,d-\)(g'B;T) £*(t;s) 



(3.42) 



In the limit R — > oo, the second and third lines are exponentially suppressed, except for the 
contribution of the massless sector of T^-\,d-\) an d the zero-mode part of E*(t;s) to the 
second line, which may be replaced by 




C* (2S) T; 



s+ 



d-1 



r(2s-l)r : 



(3.43) 
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The renormalised integral 1Z* (r^y, s) may therefore be written as 



2C(2s) lim 

T— S-oo 



T fl/2 S+! M / 



4/2 



■2£*(2s-l) lim 

T— >oo 



.7-1/2 

7-1/2 



d+\_ 



E * 



e T 2 



,1/2 



(3.44) 



where the ellipses denote exponentially suppressed corrections in the limit R — > oo. A Pois- 
son resummation over m shows that the term in round brackets is exponentially suppressed 
as T2 — > oo, and the limit T — > oo is therefore finite, and given by 



n*(r m ; S ) =RK k {r [d _ l4 _ l) -s) 

+2 £*(2s) T(2s + d-2) R 2s+d - 2 + 2 T(2s - 1) T(2s - d + 1) R^ 2s 



(3.45) 



up to exponentially suppressed corrections as R — > oo. This formula is manifestly invariant 
under ml-s, and provides the constant term for the constrained Epstein zeta series ( 3.16| 
with respect to the parabolic subgroup with Levi component SO(d — l,d — l) x R + inside 
SO(d,d). It is easy to check that it is satisfied for d = 1,2 using the explicit results in the 
previous subsection, and that it agrees with Eq. (D.16) in [20 1 for d = 3. One may also check 
that ( 3.44| l is consistent with the Laplace equation 0.8) using Eq. (A.30) in 1 15 1. 



3.4 Another modular invariant regulator 

In the context of threshold corrections in four-dimensional heterotic vacua, Kiritsis and 
Kounnas [18| have proposed a different modular invariant regularisation for on shell in- 
frared divergences, based on replacing Minkowski space by a superconformal field theory 
with the same central charge, but depending on an infrared cut-off A. For non-zero value of 
A, the threshold correction is given by a modular integral 



(3.46) 



where Z(A, t) is related to the partition function r^i u(R,r) of a compact boson of radius 
R = 1/Aby 



Z(A,t) = D • W1/A,T) , [D ■ f] (A) = 2A 2 9 A \f{2A) — /(A)] . 



(3.47) 



The integral is manifestly finite, since the cut-off function Z(A, t) decays exponentially at 
T2 — > oo. For fixed X, the integrand in ( |3.46| l agrees with the usual one in the limit A — > 0, 
since Z(A — > 0, r) = 1 up to exponential corrections. 
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To relate this prescription to ours, we apply the Ranking-Selberg method for functions of 
rapid decay, and write 



lf K {g,B) = 2Res s=1 n*(r m Z(A,r) ,s) , 



where 



K*(r m Z(A,T) ,s) = R.N. J^r {d4) (g / B;T)Z(A / T)E*(s / T) . 



(3.48) 



(3.49) 



Since Z (A, t) is of rapid decay, the sign R.N. is superfluous, however introducing it allows 
us to take the operator D out of the integral and obtain 



K k {r m Z{A,r),s) = D- R.N. J^r m {g f B;r) r (w) (l/A) E*(s,t) 



(3.50) 



Using the results of the previous subsection (with d, R replaced by d + 1, 1 / A), we find that 
in the limit where the infrared cut-off A is removed, 



n*(r m z(A,r) ,s) n*(r m ;s) 

+ 2 C(2s) £*(2s + d - 1) D • A 1 "^ 25 
+ 2£*(2s - 1) r(2s - d) D ■ A ls - d - 1 



(3.51) 



Extracting the residue at s = 1, we find that for d = 0,1, (3.461 is finite in the limit A — > 



and reduces to the standard results. For d > 2, ( 3.46 1 agrees with ( |3.20[ ) after subtracting the 
order 0(A 2 ~ rf ) divergent term. Finally, for d = 2 ( |3.46| l diverges logarithmically as A — > 0, 
and agrees with ( |3.36| l after subtracting — 21og(2eA). 



4 Modular integrals with unphysical tachyons 

In heterotic string vacua with a 'spectator' d-dimensional torus T d , an interesting class of 
couplings in the low energy effective action is determined by the one-loop integral of the 
elliptic genus |2|, 

j^r {d+k4) (g,B / y)0(r)d } i ) (4.1) 
where 

V M )(g,B,y) =t{' 1 £ ^i^! (4.2) 

PL,fR 

is the partition function of the Narain lattice, parameterised by the metric gij, two-form B« 
and by Wilson lines y\ (a = 1 . . . k), and <P (t) is a weak holomorphic modular form of nega- 
tive weight w = — A;/2 with a simple pole at the cusp, q = e 2mT = 0. This singular behaviour 
is associated with the 'unphysical tachyons' which are ubiquitous in heterotic vacua, i.e. rel- 
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evant operators in the (0, 2)-superconformal world-sheet theory which do not respect the 
level-matching condition. Although 0(t) grows exponentially as T — > ioo, the integral ( |4.1| | 
is at most power-like divergent, under the condition that one integrates first on Ti and then 
on T2. We refer to integrals of the type ( |4.1| as 'modular integrals with unphysical tachyons'. 

For simplicity we shall restrict ourselves to points in the Narain moduli space where the 
Wilson lines yf vanish, so that the lattice partition function factorises as ru + ^\(g, B,y) = 
TuA{g,B) r^Q) where Tr^ Q \ is a holomorphic modular form of weight k/2. At the cost of 
absorbing |^] -T^o) into 0, we can therefore assume that k = and is a weak holomor- 
phic modular function (i.e. a weak holomorphic modular form of weight zero with trivial 
multiplier system) 

0(r) = £ a n q n . (4.3) 

n>— k 

For physics applications, we are only interested in having a simple pole at q = (i.e. k = 1), 
however our mathematical construction works equally well for any non-negative integer k. 

For d = 0, the integrand function is holomorphic and the integral ( |4.1| | is easily computed 
by representing the integrand as a total derivative and using Stokes' theorem (2j[3|: 



/ d^( T ) = £ f 1/2 dr 1 G 2 (r)0(r) = \ U - 24 £ a n a{-n)\ , (4.4) 

where G2 is the holomorphic quasi-modular Eisenstein series of weight two, and c(n) = 
ci (n) is the sum of divisors of n: 



,m 



G 2 (T) = 1 - 24 £ = 1-24 £ a(n) f . (4.5) 

m=l ' n=l 



For d > 0, however, the integrand function is no longer holomorphic, and cannot be writ- 
ten as a total derivative. The standard approach consists instead in unfolding the lattice 
sum |3j|8]|, at the cost of obscuring the T-duality symmetry 0(d, d,Z) of the Narain partition 
function. In this section, we shall develop techniques for computing integrals of the type 
( |4.1| | while keeping this symmetry manifest. 

Our basic strategy will be to represent the weak holomorphic modular function in ( |4.1| | 
as a Poincare series, and to unfold it rather than the lattice surrj^J While straightforward in 
principle, the main difficulty in implementing this idea is the fact that the standard Poincare 



8 Note that, in general, the limit y" — > is singular due to the appearance of extra massless gauge bosons. 
However, in the case at hand these massless states do not contribute to the IR behaviour of the the modular 
integral, since we subtract their infrared divergences. Singularities do however appear at points of symmetry 
enhancement consistent with the factorisation of the Narain lattice. 

9 This is similar in spirit to the Rankin-Selberg method, however note that the Poincare series is no longer 
an auxiliary function which reduces to a constant (or a pole with constant residue) at a special value of the 
parameter, but rather represents part of the modular function to be integrated. 



19 



series representation of weak holomorphic modular functions (see e.g. 1 31 ■ 32|) 

<p{r) = \a + \ £ «nli^ E E e ^ (l-e^ , (4.6) 

-K<n<0 K ^°°|c|<K \d\<K;(c,d)=l V / 

is not absolutely convergent, so the unfolding cannot be justified. 

This problem can be circumvented if the weak holomorphic modular form <P(t) can be 
obtained as a suitable limit of an absolutely convergent non-holomorphic Poincare series (or 
a linear combination thereof), 



E( S/ K) = E(T 1/ T 2 ;s /K ) = i £ i j,i2s e ' 27rk ^- ( 47 ) 



Here the integers a and & are a solution oi ad — be = 1. The sum over c, d is absolutely 
convergent for 3R(s) > 1, and defines an automorphic form of weight 0. It may be ana- 
lytically continued to s = 0, where the summand of ( |4.7| | becomes holomorphic, and the 
Poincare series ( |4.7[ ) formally defines a weak holomorphic modular function which behaves 
as l/q K + 0(1) at q = 0. On the other hand, for any modular function E{x\, T2), the integral 
over the cut-off fundamental domain 

R t (F,s,k)= [ dpf(Ti,T 2 ) E(s,k) (4.8) 

can be computed for K(s) > 1 by unfolding the sum over c, d, and analytically continuing to 
s = as well. Choosing ¥{r\, T2) = £7^) (g, B; t), and assuming that <P(t) can be represented 
as a linear combination 

<£(T)=lim V a n E(s,-n) (4.9) 



of the Poincare series (4.12 1 in the limit s — > 0, in principle gives a way to compute the mod- 
ular integral (4.1 1 of interest, while keeping T-duality manifest. To carry out this program, 
we shall compute the Fourier expansion of ( |4.7| , and study its limit as s — > 0. 

Before doing so however, let us first note that the non-holomorphic Poincare series ( |4.7| | 



is not an eigenmode of the hyperbolic Laplacian 

A = 2 if 9 t 9 T , (4.10) 
but, rather, satisfies [33 1 

[A + is(l-s)] E(s,k) = 2tcksE{s + \,k). (4.11) 

This equation in principle allows one to determine the analytic properties of E(s,k) for 
3f£(s) < 1 from the knowledge of the spectrum of the Laplacian A. 
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4.1 Fourier expansion of the non-holomorphic Eisenstein series 



The Fourier expansion of E(s, k) can be computed by standard methods |33 -37J. After ex- 
tracting the term with c = 0, d = 1, and setting d = d' + nc in the remaining sum, a Poisson 
re-summation over n yields 



E(s,k) = t| e 



s „—1n\KT 



Y_ E n (s,K)e 



ImtiT 



(4.12) 



where the n-th Fourier coefficient 

E n (s,K) = rt s £ S(n '- K;c) 



~2s 



A n {T lr c;s,k) 



(4.13) 



c=l 



is expressed through the integral 

f+ ca 

A n (r 2 ,c;s,K) = d£(r + l)~ s exp 

J — CO 

and the Kloosterman sum 



27Ti 



c 2 r 2 (t 



nr 2 (t + i) 



S(a,b;c) = ^ exp ( ^^-(a d + bd 1 ) \ , 

de(Z/cZ)* x L 



(4-14) 



(4.15) 



where d 1 stands for the inverse of d modulo c (see AppendixjXJfor some relevant properties 
of the Kloosterman sums). The Kloosterman-Selberg zeta function 



^ S(a,b;c) 



(4.16) 



oo 



which is holomorphic for 3ft (s) > 1 and admits a meromorphic continuation to all complex 
values of s, will play a central role in what follows. 

The integral ( |4.14| | can be computed by Taylor expanding exp ^ &r™+i) ) 1 anc ^ usrn § * ne 
standard formula for the integral 



/oo 
(t + i)-« (t-i)-Pe- 2niut dt = < 
-CO 



22-"-^ r(tt+/3-l) 



if w = 0, 



(4-17) 



valid for u G R, 3ft(a + /3) > 1. Here, the function cr{rj,(x.,fi) is related to the Whittaker 
function W^ /V (x) by 



'fa, cc, p) = r(fi) y-^W 2 e>1 /2 aJ± ± fa) . 



(4.18) 
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In this fashion, we find that the constant term is given by 



E (s,k) = £ 
m=0 



2 2(l- S ) 7r(7rK) m r(2s + m _ 1) 

— Z(0, — k; s + m) t 2 



m! r(s + m) r(s) 
while the positive-frequency Fourier coefficients are given by 



(4.19) 



00 (27i) 2 ( s+m ) n 2s+m ~ l K m 
E n>0 (s,K) = V ± — , w . w — cr(47rnT 2/ s + m / s) Z{yi,-k)s)t{ (4.20) 
£± m\r(s)r(s + m) 



and the negative-frequency Fourier coefficients (besides the first term in (4.12 1) are given by 



~ (27r) 2 ( s+m ) (- n ) 2s+m - l K m 
~ m\r{s)r{s + m) 



g inm 2 a (—4777^ S/ s _|_ m ) Z(tl, —K)S + m) t| . 

(4.21) 



Using standard estimates on the Whittaker function and on the Kloosterman-Selberg zeta 
function, one may check that all these series are absolutely convergent for 9ft (s) > 1, and 
therefore E(s, k) is a meromorphic function of s, analytic in the half-plane 3ft (s) > 1. 

It will be useful to rewrite these expressions in a different form. First, using the fact that 
the Kloosterman-Selberg zeta function Z(a,b;s) at a = can be expressed in terms of the 
Riemann Zeta function via ( A.4| , the constant term ( 4.19| l may be written as 



Eo(s,k) = £ 



~ 2 2 (!- S ) 71 (7IK) m r(2 S + m-l) <7l-2s- 2m (K) 1 



.„_ m\r(s + m)r(s)Z(2s + 2m) 



(4.22) 



It may be checked that (422) satisfies the differential equation ( 4.11 1 with E(s,k) replaced 
by Eo(s,k). Second, using the standard series representation of the Whittaker and Bessel 
functions 



I (z) - f {Z/I)2m+V 
vK) ^ m\r(m + v + l)' 



(4.23) 



one may carry out the sum over m and obtain, for n > 0, 



E n ( S/K )-2zr^ r(s)p , ^ 



E 

c>0 



S(n, -k;c) 

c s+l+p 



^s— p— 1 



4tt ,—\ 

C J 



(4.24) 



As an example, we see that the Poincare series E(s,k) can be analytically continued to s = 
1/2, where its Fourier expansion becomes 

E(J,k) = v^ e - 2 ^ + \2^t 2 a (K) + ^y^^ + ^y^JV + ■ ■ ■ J + ■ ■ ■ ■ < 4 - 25 ) 



For K = 1, this reproduces Eq.(3.45) in (38j. 
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4.2 Holomorphic limit 

Let us now study the limit s — > of the non-holomorphic Poincare series E(s,k). At this 
value, the summand in ( |4.7| becomes a holomorphic function of T. Thus, one expects that 
the analytic continuation £(0, k) will be a holomorphic modular function^ In this limit, the 
constant Fourier coefficient ( |4.22| | reduces to 

E (0, k) = 12 £t(jc). (4.26) 



Moreover, the negative frequency Fourier coefficients (4.21 1 (besides the first term in ( 4.12| ) 



vanish in this limit, while the positive frequency Fourier coefficients become 

« / \ ~ / K «-i S(n. — k;c) t /47T , — \ .„ „_ 

E„(k) = 2tt J- £ ' j Ii — • (4.27) 



Thus, the analytic continuation of E(s, k) at s = is given by 



E(0 /K ) = <f * + 12cr(K) +2tt E \P E SK K;C) Ii f — v 7 ^ ) ^ • (4.28) 



For k = 1, this is recognised as the Petersson-Rademacher formula 1 39 40 1 for the Klein 



modular invariant /(t) = 1/q + 196884^ + . . . , up to a suitable additive constant, 

E(0,1) =;(t) + 12. (4.29) 



For k > 1, the r.h.s. of (428) is obtained from its value at K = 1 by acting with the the Hecke 
operator 



(t,.#)(t)= E E *(°^)- 

«,d>0 fcmodi V " / 



(4.30) 



which acts in the space of weak holomorphic modular functions and maps = 1/q + 0(q) 
to T K ■ = l/q K + 0(q). This can be checked using Selberg's identity ( |A.2| for the Klooster- 
man sum and the fact that T K acts on the Fourier coefficients as 

(T K -E) n = K E d~ l E nK/d 2. (4.31) 
d\(n,K) 

Thus, the analytic continuation of E(s, k) at s = is, in general, the weak holomorphic 
modular function 

E(0,k) = T K -;(T) + 12cr( K ). (4.32) 



10 For negative modular weight, analytic continuation can lead to holomorphic anomalies but, as we shall see, 
these do not occur for w = 0. We thank J. Manschot for discussions on this issue. 
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Indeed, a numerical evaluation of ( 4.28[ ) reproduces the known coefficients of the /-function 
and its images under T K . 



4.3 Shifted constrained Epstein zeta series 

In this subsection we shall evaluate the modular integral ( |4.8| l for an arbitrary modular func- 
tion F(ti, t 2 ) of moderate growth. In the interest of simplicity, and with a view towards our 
main goal ( |4.1| , we assume that 



Fo(t 2 ) = t| /2 , (4.33) 
and we require that there exists an a > such that 

F K (T 2 )e 2nKT2 ~ e- aT2 (4.34) 

as t 2 — >• 00. Here F k (t 2 ) = X^i/ 2 f( T l/ T 2) e _27nKT:1 dii is the K-th Fourier coefficient of F. 
Notice that these conditions are satisfied by the lattice partition function Fud) away from 
points of extended gauge symmetry. 

By the same reasoning as in eqs. ( |2.18| -( 2T20) , for large enough 3?(s) we can unfold the 
sum over (c,d) in the non-holomorphic Poincare series ( |4.7| |, arriving at 

R t (F,s,k) = f dT 2 r 2 s -V^F K (T 2 )- / dpF(r lf r 2 ) (£(s, K )-T| e - 2 - KT ) . (4.35) 

JO J J-—J--J- 

By our assumption ( |4.33[ ), the term on the first line is finite as T — > oo. The term on the 
second line can be decomposed into 



- / dp (F^T^-if 2 ) (E(s,K)-r| e - 27riKT )- £ e m (s, K ) / dr^ 1 s m+2 , 

(4.36) 

where 



p ( Q r \ _ ? 2(l-s) - ( nK ) m <Tl-2s-2m{K) F(2s + ffl - 1) 

m{ ' j m!F(s + m)F(s)C(2s + 2m) " 



(4.37) 



The first term in this expression is exponentially suppressed in the limit T — > oo, thus, up to 
exponentially suppressed terms, the modular integral ( |4.8| l is given by 

R t (F,s,k)= dT 2 T|-V^F K (T 2 )+£^ e !^ K l + .... (4.38) 

Jo - ~ \-s-m 

For 3?(s) > d/2, the second term in ( |4.38| l is suppressed as T — » oo, and the cut-off integral 
( |4.8| l converges to the first term in ( 4.38| l as the cut-off is removed. We define the renormalised 
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integral as the analytic continuation of this result to all s, namely by the Mellin transform 



R.N. J F(zi, t 2 ) E(s, k) dp = dr 2 T 2 s - 2 e 27reT2 F k (t 2 ) 



(4.39) 



As in the case considered in Section 2, the analytic structure of (4.39 1 can be read off from the 
growth of e 2nKT2 F k (t 2 ) near t 2 — > oo. Unlike the case studied in Section 2, the integral (439) 
is not expected to obey a simple functional equation under s — > 1 — s. 

Let us now specialise to the case F(ti, t 2 ) = F(^v The K-th Fourier coefficient is given 

by 



F K (r) = r 2 d/2 £ e - 7rT2 ^)^(p 2 L -p|-2 K ) / 



(4.40) 



Pl,Pr 



where the sum is restricted to lattice vectors satisfying the usual level-matching constraint 
Pl~ Pr = 2k. For such vectors, v\ J r'p\ = 2 {Pr + jc), and the condition ( 4. 34) is therefore 
obeyed away from loci on the symmetric space 0(d, d, R) /O(d) x O(d) where p\ vanishes 
for some lattice vector satisfying the constraint above. Physically these loci correspond to 
points of enhanced gauge symmetry where additional massless states occur, leading to new 
infrared divergences. We shall always work away from such points. 



Under this assumption, the Mellin transform (4.391 can be computed explicitly by inte- 
grating term by term, as in ( 3.14} . The result is a 'shifted constrained Epstein zeta series' 



R.N.J^r {did) (g,B)E(s, K ) 



r(s + 1 -l) 

TT'2 1 



where 



£ v {g,B;s + \d-l,K), 



(4.41) 



S(pl-pl-2K) 



^pApi + pI- 2 *) 



Pl,Pr 



(4.42) 



This series is absolutely convergent when K(s) > d. Unlike the 'unshifted' Epstein zeta 
series ( |3.15} , which we henceforth denote by £y (g, B; s, 0), the shifted series ( 4.42} diverges as 



(2p| ) s at the points of enhanced gauge symmetry where the norm y\ of some lattice vector 



satisfying the constraint vanishes. Moreover, using the differential equation ( 3.10a} , satisfied 
by the lattice partition function, together with the differential equation ( 4.11} satisfied by the 
non-holomorphic Poincare series, we find that the shifted constrained Epstein zeta series 
satisfies 



A 



SO(d,d) 



s(s-l) + J <*(<*- 2)1 £i{g,B;s + \d-l,K) = 2TCKs£${g,B;s + \d,K) . (4.43) 



In the limit s — > 0, the r.h.s. vanishes and therefore £y (s + ^ — 1, k) is an eigenmode of the 
Laplacian on 0{d, d, R)/0(d) x O(d) with eigenvalue d(2 — d)/A. This argument assumes 
i, k) is finite at s = 0, otherwise the r.h.s. of ( |4.43} may be non-vanishing. 
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4.4 Holomorphic limit of the integral 

Let us now investigate the limit of the integral ( |4.41| | at s = 0. As discussed below eq. ( |4.7| |, 
the non-holomorphic Poincare series becomes holomorphic in this limit, and one may hope 
to recover a modular integral of the form ( |4.1[ ). Of course, infrared divergences require spe- 
cial care. 

Returning to (438) at finite infrared cut-off T, we see that all power-like terms vanish 
in the limit s — > 0, appart from those corresponding to m = 1 (whose coefficient goes to a 
constant e\(s,K) = 12a (k), c.f. ( 4.29| ) and torn = d/2 (when d is even), whose coefficient e m 
vanishes linearly in s but which are divided by a vanishing number. Considering first the 
case where d is odd and greater than 3, we thus have 



lim 



d ,/ £ ( s, k) f ( „.,, , ) = r(z ,_ i 1) £y{g,B;\d — 1,k) + 12a( K ) T ' 



-l 



d _1 
2 1 



(4.44) 



s->0 / 

On the other hand, in the limit s — > 0, E(s,k) reduces to the holomorphic modular forms 
T K ■ ;'(t) + 12 ct(k), cf. ( |4.32[ ). The constant term in this expression is responsible for the 
divergent term in ( |4.44| . After subtracting this divergence, we conclude that 



RN.j^dp [T k .;(t) + 12(7-00] r (w) 



Q|^4(g,B;§d-l,K), 



71" 



(4.45) 



or, using ( |3.14| l, 



R.N. / dfiT K -j(T) r m 



r(i-i) 



71 



d/2-1 



£ a v {g, B;\d — 1,k) — 12 cr( K ) (g, B;\d — 1, 0) 



(4.46) 



If d is even and d > 4, there is an additional constant term on the r.h.s. of ( 4.44| |, coming from 
m = d/2: 



I 
2 

After subtracting the divergent term, we conclude that for d > 4 even, 



55 (...) - -^rzr ft 1 <« - W + in + • (447) 



R.N 



t) r, 



(d,d) 



r(f-i) 



7T 



•d/2-1 



B; i d - 1, k) - 12 o-(k) 5# (g, B; \ d - 1, 0) 



87T(7TK) d/2 (7- 1 _ (? (K) 

(d/2)!(d-2K(d) 



(4.48) 



In the remainder of this subsection we deal with the low dimensional cases. 
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For d = 0, both m = and m = 1 contribute, leading to 

I dp [T K • ;(t) + 12<t(k)] = -12ct(k)T _1 - 4tt^(k) + . . . . (4.49) 
Using fj. dp = j — Y' we arrive at 

/ dp T K • ;(t) = -8ttc7(k) , (4.50) 

which agrees with the result of Stokes' theorem ( |4.4| . Notice that this result is independent 
of the cutoff T, since the only contribution to the modular integral comes from the region 
t 2 < 1 of T. 

For d = 1, only m = 1 contributes, leading to 

/ dp [T k -;(t) + 120-(jc)] r (u) (R) = -24a( K )T- 1/2 - An £ |m/R-nR| (4.51) 

3~T m,n>0 

mn=K 

and therefore 

/ dpT K .;'(T)r (11) (R) = -4zr £ |m/R-nR| - 47rcr( K )(R + 1/R) . (4.52) 



m,n>0 
mn=K 



For R > a/k, Eq. ( 4.52| | reduces to — 87rc(jc)R, which would be the result of naively applying 



the unfolding trick to the partition sum ^(R), as in |9J. Similarly, for R < 1/yK, Eq. 



( 4.52 | reduces to — 8tcct(k) /R, as required by T-duality. For 1/ y/x < R < yjx, however, the 



naive unfolding of (R) fails, while our method still applies. We note that the Eq. ( 4.52| | 



for k > 1 can be derived from the k = 1 result by observing that the lattice partition function 
satisfied"! 

T K -r {u) {R)r) = y/x £ r (w) (RVm7n;T), (4.53) 

m,«>0 
mn—K 

and invoking the Hermiticity of the Hecke operator T K with respect to the Petersson product. 
Eq. ( |4.53| maybe viewed as a p-adic analogue (for k = p prime) of ( 3.10a[ >. 



Finally, we consider the most complicated case d = 2. In this case the coefficient e\ 
behaves as 

ei = a + /3s + 0(s 2 ), (4.54) 

with 

« = 12£t(k) / /S = -24 ct(k) [7- 121ogA + log(47r)] -24jcc/ 1 (k) , (4.55) 



11 This identity is easily established for k a prime number, and can be extended to the general case using the 
Hecke algebra T K T K i = Y,d\(K k') d T KK * /gi (note the non-standard normalization of the Hecke operators in l |4.30| l). 
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where A is the Glaisher constant. After integrating over x 2 , we arrive at 

/ dfiE(s /K )r i2i2) = T ^-£ v (g,B;s, K )-- + cclogT-B + 0{s). (4.56) 

J J- j- 7T S 

The pole and logarithmic divergences on the r.h.s. originate from the constant term 12 c{k) t| 
in the limit s —¥ of E{s, k). After subtracting these terms, we conclude that 



/ dy, T K ■ j(r) r (2/2 ) =lim 



£y{g,B;s,K) 



B - 12 a( K )I 2 . (4.57) 



7T S 

In particular, we conclude the Epstein zeta series £ v (g, B; s, k) has a zero at s = 0. 

On the other hand, the integral ( 4.57| may be determined using harmonicity and the 



singularity structure near points of enhanced symmetry. E.g. for K = 1, states with momenta 
mi = tlx = 0, ni2 = —n-2 = ±1 become massless^Jat T = U, and induce a singularity of the 
form 

-41og|r-tf|. (4.58) 

More generally, massless states arise whenever j(T) = j(U), as required by invariance under 
the T-duality group SL(2,Z)7 x SL(2,Z)u k Combining this singularity structure with 
harmonicity, one can uniquely determine 

d ¥ ;(t) T m (T, U) = - log \j(T) - j(U) | 4 + const . (4.59) 

This result was confirmed in |8) by using the unfolding trick and of Borcherds' product 
formula for j(T) — j(U). 

Similarly to the one-dimensional case, the integral ( 4.57| | can be evaluated for k > 1 by 



observing that the lattice partition function satisfies the p-adic analogue of ( |3.39| l 

Tk^ r (2,2) = T K T T( 2/ 2) = T( 2/2 ) , (4.60) 



where T^' U ^ are the analogues of the Hecke operator ( 4.30| l, now acting on the moduli T and 



U. Invoking the Hermiticity of T K with respect to the Petersson product, one finds that ( 4.59| 
generalizes into 

J^d ¥ T K -j{r)r m {T,U) = -i(Ti T) + Ti U) )log \j(T) -j(U)\ 4 + const. (4.61) 

This result is consistent with the fact that additional massless states arise at T = k U, together 
with the images of this locus under T-duality. 

Comparing ( 4.61| with ( 4.57[ ), we arrive at interesting identities between the constrained 



Epstein zeta series £y(g,B;s,K.) at s = and the function. It would be interesting to com- 



12 



Note that 4 additional massless states appear at the special values T = U = i and T = U = p, where 



p = e in/3 . 
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pute the constrained Epstein zeta series directly, perhaps along the lines of Section 3.2 and 
derive these identities independently. 

More generally, it will be interesting to extend our method to modular integrals of non- 
symmetric lattice partition functions F^_|_jt /( j. While the generalisation of the non-holomorphic 
Poincare series ( |4.7| to negative weight is obvious, its analytic continuation to the relevant 
value of s where the summand becomes holomorphic is subtle, and involves an interesting 
interplay of holomorphic and modular anomalies (37|. We hope to discuss this problem in 
future work. 
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A Properties of Kloosterman sums 



Kloosterman sums play an central role in Number Theory. The classical Kloosterman sums 
for the modular group T = SL(2,Z) are defined by 



S{a,h;c) = exp 

de(z/cZ)* 



2 m 



(ad + bd- L ) 



(AT) 



where a, b and c are integers, and d^ 1 is the inverse of d mod c, and enter in the explicit 
expression of the Fourier coefficients of modular forms. S (a, b;c) is clearly symmetric under 
the exchange of a and b. Less evidently, it satisfies the Selberg identity 



S(a,b;c)= £ dS(ab/d 2 ,l;c/d) 

d\gcd(a,b,c) 



(A.2) 



In the special case a ^ 0, b = 0, the Kloosterman sum reduces to the Ramanujan sum 

S{a,0;c) = S(0,a;c) = £ exp ( — ad] = £ dji(c/d), (A3) 

d€(Z/cZ)* V C ' A\ a rMra\ 



d|gcd(c,n) 



with }i(n) the Mobius function. For a = b = 0, S(a,b; c) reduces instead to the Euler totient 
function (p(c). 
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We now turn to the Kloosterman-Selberg zeta function ( 4.16| . Using the trivial bound 



\S(a, b;c) \ < c, one sees immediately that the sum over c converges absolutely when 3?(s) > 
1. The Weil bound \S(a,b;c)\ < 2 V (°> ^Jcgcd{a,b,c), where v(n) is the number of divisors 
of n, shows that Z(a,b;s) is in fact analytic when K(s) > 3/4. When one or both of the 
arguments vanish, it can be expressed in terms of the Riemann zeta function via 

z(a0;s) = ^T' z(a±K;s) = ^f {K ^ 0) ' (A ' 4) 

where c s (n) is the divisor function 

a s {n)=Y J d s . (A.5) 

d\n 

The usual notation cr{n) = (T\{n) is used throughout the paper. 
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